Spin Content of A in QCD Sum Rules 
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We calculate the isoscalar axial-vector coupling constants of the A hyperon using the method of 
QCD sum rules. A determination of these coupling constants reveals the individual contributions of 
the u, d and the s quarks to the spin content of A. Our results for the light-quark contributions are 
in agreement with those from experiment assuming flavor SU(3). We also find that the flavor-SU(3)- 
breaking effects are small and the contributions from the u and the d quarks to the A polarization 
are negatively polarized as in the flavor-SU(3) limit. 
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I. INTRODUCTION 



According to data from polarized lepton-nucleon deep 
inelastic scattering (DIS) which was initially reported by 
the EMC Collaboration [l| and confirmed subsequently 
by several other experiments 0, 13, 0] , only a small frac- 
tion of the nucleon spin is carried by the valence quarks. 
This observation has triggered many research activities 
and puzzles in understanding the spin content of the nu- 
cleon (see e.g. Refs. 0, [1| for a review). In this respect, 
it has been realized that the polarized A, having all spin 
carried by the s-quark while the u-d quark pair is coupled 
to 5' = 0, / = 0, provides a special example in the naive 
quark model. Contrary to the naive expectations, the in- 
terpretation of the experimental data together with the 
SU(3) F symmetry is that ~ 60% of the A spin is carried 
by the s (and s) quark, while ~ —40% originates from 
u (and u) and d (and d) quarks p]. One important as- 
pect of this interpretation is that the u-d quark pair has 
negative polarization. Actually, this result is supported 
by some model-dependent approaches as well (see Ap- 
pendix An investigation of this interesting issue has 
the potential to shed light on the "spin crisis" and there- 
fore attracted considerable attention 0, H, [1, [13) E!, ■ 
Experimentally, the polarization of A is of special interest 
because it can be easily measured from the nonleptonic 
decay A^PTT [H, [11, [M [13 ■ 

One open question in this framework is how sensitive 
the A spin structure to SU(3)i? breaking effects is. While 
it is claimed that the SU(3)i? breaking may lead to a 
change in the sign of the u- and the d-quark polariza- 
tions in A lattice [l^ and some model-dependent 
works find that A is insensitive to SU(3)i^ break- 
ing. The isoscalar {g\ and g\), octet {g1) and singlet 
((^q) axial- vector coupling constants of A (shown gener- 



ically as gA = gA{q^ = 0) throughout the text) can be 
expressed in terms of the fractional contributions of the 
quark flavors, Aq, to the A spin content as 



Ad, 



g\ = Au + Ad-2As, 



9l - As, (1) 
gj^ = Au + Ad + As, 



in the SU(3)f limit. Therefore, a determination of these 
coupling constants reveals the spin content of A. 

Our primary aim in this paper is to calculate the 
isoscalar coupling constants f/^ and g% . For this purpose 
we use the method of QCD sum rules (QCDSR) [H [H, 
[20l . [2l|. Note that this is reminiscent of the works in 
Refs. m [13, S [H [13, m, where the axial- vector 
coupling constants of the nucleon have been calculated. 
Using this method, we also extract the contributions of 
the u and the d quarks to A spin content in the SU(3)_f- 
breaking case in order to see how sensitive the results are 
to symmetry-breaking effects. 

We have organized our paper as follows: In section HIl 
we derive the QCD sum rules for gA and give numer- 
ical analysis and the results in section IIIII Finally we 
conclude in section ITVl 



II. FORMULATION 

We start with the correlation function of two A inter- 
polating fields in the presence of an external constant 
axial- vector field Z^, defined by 



ii{p)-^z-n^{p) + o{z'). 



(2) 



This correlation function is computed by adding the term 
A/: = -^3,9^759, (3) 
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to the usual QCD Lagrangian, where gq is the coupling 
of the quark field to the external field and we use the 
notation ^= Z^^^. 
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The most general A interpolating field is defined as a 
mixture of two independent local operators via the mix- 
ing parameter t: 



VA = 



ea6c|2 [ulCj5db] Sc + [ulC"f5Sb] dc 

- [dlCj^Sb] u, + t{2 [u^Cdb] j5Sc 



754 - [dlCsb] 75Uc) |, 



where a, b, c are the color indices, T denotes transposi- 
tion and C = 47^7". The choice t = —1 gives the loffe's 
current, which is often used in QCDSR calculations. In 
our numerical analysis, we take t = —1.2 which produces 
the optimal interpolating field [29]. In Eq. Il{p) is 
the correlation function when the external field is absent 
and corresponds to the function that is used to determine 
the A mass, while 11^ {p) represents the linear response 
of the correlator to a small external axial- vector field Z^. 
In the presence of an external axial-vector field, Lorentz 
invariance of the vacuum is broken and new vacuum con- 
densates appear as 



{qgcG^ul^q) z ^ 9qZ^K (qq) , G^i. 



1 



(6) 



which are defined in terms of the susceptibilities x ^-nd n 
with the QCD coupling-constant squared — Airag- 

We can bring the correlation function in Eq. ^ into 
the form 



n^(p) = nf(p2)zz^a''>,75 

+ n^{p^)z-p^^5 + Til{p^) 



(7) 

(8) 



The Operator Product Expansion (OPE) sides of the sum 
rules are obtained by inserting the interpolating field (j4|) 
into the correlation function ^ and evaluating the time- 
ordered contractions of quark fields, which include the 
quark propagators [13, [H, [2^. On the OPE side, we 
include the terms up to dimension 8. The perturbative- 
correction terms in order of as may become important 
especially at lower Borel-mass region but they are ex- 
pected to give smaller contribution to the sum rule we 
choose to work with (see Section IIIip '2P] . Therefore 
these correction terms are neglected in this work. The 
phenomenological side is obtained via a dispersion re- 
lation, which is written in terms of hadron degrees of 
freedom. Finally, the QCD sum rules are constructed 
by matching the OPE sides with the phenomenological 
sides and applying the Borel transformation. The details 
of this procedure can be found in the extensive litera- 
ture on QCDSR. Omitting the details, here we give the 
final forms of the sum rules for gA as obtained at three 



different Lorentz-Dirac structures: 

Ci ag M'^ Eo L^l'^ + C2msX aq Eq L'^^^ 

+ C5Xml al L-^^l'^'^ -f- Q a, & -|- G7 a 
X 



(9) 



A* mA 



z -pit^b ■ 

C[ A/6 El + c'2xag Eq L'^/^ 



(5) ^75 



+ C!,xagbL-^/^ + C;,m,xa'gL-^^' 
+ C'm 1^"^/"'' 
: (gA + A'M^), 



(10) 



Ci M« E2 L-^^^ + C2Xag E^ L'^'"" 

■^C^nag Eo i-68/8i + q b Eq L'^/g 
+ C5 rus ag M'^ Eo L'^^^ + Cq L'^/^ 
+ C7xagb L-^l^ -f Cg m, X al L^^/^ 



(11) 



= (5A+iM2), 

where we have defined 

[t-l) 



XI (M2 - 2ml) 



Ci 
C2 

C5 
C7 



18 

12 

(t-i) 



18 

+ 2.g,(t + 5)/ 
{t-l) 



igu+gd){!it + l)+gs{t + 3)f 
{9u+9d), 

{9u+9d)[{f-5)t-{f + l)] 



(g„-f.gd)[(13/-47)t+(5-13/)] 



324 

+ 2g447~5t)f 

it -I) 
144 

+ 2g,(t + 5)/ 



(12) 



{9u+9d)h{3f + 7)t-i5f + 3)] 



(t-l) 



1296 

(i-1), 



igu+9d)int + 7)+9si35t + 37)f 



54 



i9u+9d) 2{t-l) + i5t + l)f 



3 



r' - — 



Ci = 

C2 = 



{hf + 26i + 5)(.g„ + gd) + (llt^ + Ut + ll)^, 
- 2(<2 + 7f + l)(g„ + gd) + {t- Iffgs 
(2M2 + 50i + 29)(5„+5d) 



1 

"T08 

+ 2(13i2 + 10i+13)/5, 
1 r 



96 



[t - lf{gu + gd) + (13<2 + lot + 13)5, 
1 r. 



(13) 



1 

432 
1 

9 
1 

243 



(5i + l)[4(t - 1) + 3(t + 5)/](g„ + .grf) 

(10/t + t + 2/-l)(.9„+gd) 
+ (2/f + 33t+ 10/ + 39)5,1, 

(7^2 + 4t + 7)(5„ + 5d) + (t- 1)^5, 



36 
(t-1) 
27 



(i" + 7< + l)/(5,+5d) 



[(7/ + 10)t2 + (13/ - 8)f + 7/ - 2](5„ + gd) 



1 

24 
1 

18 



^ 36 



6*4 = — 

96 



(bf + 26t + 5)(5„ + 5d) + (IK^ + 14i + 11)5, 
- (5i2 + 8i + 5)(5„+5d) 
2(19t2 + 16t+ 19)/5, 
[(25t2- 14t + 25](5„ + 5d) 
2(41<2 + 26t + 41)/5, 
{t ~ 1)2 (5„ + gd) + {lif + m + 13)5, 



(5t + l) 

12 
(t-1) 



27 



[5/ + <(/ + 4)-4](5„+5d) 
[(4/ + t(20/-l) + l](5„+5<i) 



[i(33 - 2/) - 10/ + 39]5, 



(14) 



1 



{bf - t + h\{gu + gd) 



216 

+ 2(10*2 + 7i + 10]/5. 



1 

18 
1 

162 



[(3/ + 10)i2 + 4(3/ - 2)t + 3/ - 2](5„ + gd) 
[(7/ + 10)*^ + (13/ - 8)t + 7/ - 2](5„ + gd) 



Here Af is the Borel mass and the overlap ampUtude is 
defined via (0|?7a|A(p)) = Aaw(p) [^{p) is the Dirac spinor 
for A with momentum p] with = 327r^AA- We have 
also defined the quark condensate ~ ~ {2Tr)'^ (qq) , and 
the quark-gluon-mixed condensate (g5cf • Gq) ~ rn^ (qq) . 
The flavor-symmetry breaking is accounted for by the 
factor / = (ss) / (qq) . The continuum contributions are 



included via the factors 

En = l- (l + x^ 



(15) 



with X — vP'/M'^, where w is the continuum thresh- 
old. The corrections that come from the anomalous di- 
mensions of various operators are included with the fac- 
tors L = log(M2/AQ(^^)/ log(^2/j\^^^)^ where ^ is the 

renormaUzation scale and Aqcd is the QCD scale pa- 
rameter. 



III. RESULTS 

In principle one can use any of the three sum rules 
to calculate gA, however, not all of them work equally 
well due to continuum effects and insufficient OPE con- 
vergence. For the calculation of the nucleon axial-vector 
coupling constants, the sum rule at the structure Z ■ p^ 
75 has been favored over the others in the literature. On 
the other hand, it has been found in Ref. ^] that this 
sum rule fails to have a valid Borel region whereas the 
sum rule at the structure iZ^ cr^'^ Pi^j5 satisfies OPE con- 
vergence and pole dominance, therefore has a valid Borel 
window. 

The valid Borel regions are determined so that the 
highest-dimensional operator contributes no more than 
about 10% to the OPE side, which gives the lower 
limit and ensures OPE convergence. The upper limit 
is determined using a restrictive criterion such that the 
continuum-plus-excited-state contributions are less than 
about 30% of the phenomenological side, which is im- 
posed so as to warrant the pole dominance. Using these 
criteria we have similarly found that the sum rule in 
^ has a valid Borel window while those in (fTO|) and 
(fTTj) are seriously contaminated by continuum contribu- 
tions. Therefore we choose to work with the sum rule 
at the structure iZ^a^^ p^^^ in ([9]). In order to ob- 
tain the corresponding sum rules for 5^ {9 a) we set 

gu = gd = ^ {gu = gd = 0) and 5, = (5, = 1). 

We determine the uncertainties in the extracted pa- 
rameters via the Monte Carlo-based analysis introduced 
in Ref. [2^. In this analysis, randomly selected, Gaus- 
sianly distributed sets are generated from the uncer- 
tainties in the QCD input parameters. Here we use 
3 b = {glG^) = 1.2 ± 0.6 GeV^ 



a, = 0.52 ± 0.05 GeV 



TO^ = 0.72 ± 0.08 GeV^, and Aqcd = 0.15 ± 0.04 GeV. 
The flavor-symmetry breaking parameter and the mass of 
the s-quark are taken as / = (ss) / (uu) = 1 and m, = 0, 
respectively, in the SU(3)f limit. In the SU(3)i?-broken 
case, we take these parameter values as / = 0.83 ± 0.05 
and m, = 0.11±0.02 GeV. The values of the vacuum sus- 
ceptibilities have been estimated in Refs. [2^ [27l. [sol. l3l| . 
We consider the values x^tg — 0.60 GeV^ and KOq = 
0.05 GeV* for 5^ and 5^. The continuum threshold 
is taken as w = 1.5 GeV in the SU(3)f limit and as 
w = 1.7 GeV in the SU(3)i?-breaking case. 
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FIG. 1: The left-hand (sohd red curves) and the fitted right- 
hand (dotted lines) sides of the sum rules Q for g\ and 
in their valid Borel regions. The bands show the errors as ob- 
tained from the Monte Carlo-based analysis and the diamonds 
mark the values based on expectations from experimental re- 
sults assuming SU(3)f. 
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For normalization of the sum rule ^ , we use the chiral- 
odd mass sum rule, which is obtained using the invariant 
function Ili{p^) as follows: 

C'l aq El + C2 ml a, Ea + C3 a, b 

+ 6*4 rris E2 + C5 ms a\ + Cq 7n, b Eq ^-^g^ 



where 
Ci 
C2 
C3 



12 

24 

{t-l 



288 
12 



11/+ 10+ (13/ + 2)i 
7/ + ll + (ll/ + 7)t 
13/ + 2+(ll/ + 10)i 
(13t + ll), 



(17) 



1 r 
18 



3(5i2 + 2i + 5) + (i-l)(i + 5)/ 
(t-l) 



96 



-(lit + 13). 



Note that the chiral-odd mass sum rule has been found 
to be more reliable than the chiral-even one [231 . 



TABLE I: The isoscalar {g\ and g\) axial-vector coupling 
constants of A as obtained from QCDSR. For comparison, we 
also give the coupling constants from naive SU(3)f assum- 
ing (A'^|s7;j75s|A/') = (denoted by SU(3)F[naive]), and those 
from DIS data assuming SU(3)f (denoted by SU(3)f[DIS]), 
which are obtained by inserting E = ?>F — D and E = 0.33 
(central value as reported by the HERMES Collaboration [3|), 
respectively, into the SU(3)f relations. As for the SU(3) f pa- 
rameters we use F/D = 0.575 and _F + D = 1.269 



9A 


SU(3)f [naive] 


SU(3)f[DIS] 


QCDSR 




-0.15 


-0.32 


-0.37 + 0.13 


Qa 


0.73 


0.65 


0.51 + 0.11 



We first concentrate on the sum rules for the isoscalar 
coupling constants g\ and in the SU(3)i? limit. The 
Monte Carlo-based analyses of the sum rules are per- 
formed by first fitting the mass sum rule (|16p to simulta- 
neously obtain m\ and Aa , and the obtained value for the 
overlap amplitude is used in the sum rules of gA for each 
corresponding parameter set. In Fig. [2 we plot the left- 
hand and the fitted right-hand sides of the sum rules ([9]) 
in their valid Borel regions. The bands show the errors 
as obtained from the Monte Carlo-based analysis. 

Our numerical results are given in Table IH For com- 
parison, we also give the couphng constants from SU(3) p 
assuming (7V|s7^75s|7V) — 0, and those from DIS data as- 
suming SU(3)f relations given as g\ — 2/3(S — D) and 
g\ — 1/3 (S + 2D), where S is equivalent to the flavor- 
singlet axial- vector coupling constant, gj^. The value we 
obtain for 5^ in QCDSR, namely g\ — —0.37 ± 0.13, is 
in nice agreement with the experimental result, while g\ 
lies slightly lower. 

FIG. 2: Same as Fig.[T]but for g\ in the SU(3)F-broken case. 
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The SU(3)F-breaking effects are accounted for by 
restoring the physical values of the parameters ms and / 
in the sum rules (O and . We consider only the sum 
rule for g\ since the susceptibilities associated with this 
coupling are unaffected with SU(3) p breaking. We apply 
a similar procedure as above where we obtain Aa from the 
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mass sum rule in (|T6)) with the SU(3)i?-breaking effects 
and this value is used in the sum rules of gA for each cor- 
responding parameter set. In Fig. [21 we plot the left-hand 
and the fitted right-hand side of the sum rule ^ for 
in the SU(3)F-brokcn case. We obtain g\ -0.29±0.22, 
which is consistent with the value obtained in the SU(3)_f 
limit. This implies that the SU(3)_F-breaking effects are 
small and the contributions from the u and the d quarks 
to the A polarization are negatively polarized as in the 
SU(3)f limit. 



IV. CONCLUSION 

In conclusion, we have calculated the isoscalar axial- 
vector coupling constants of A using the method of 
QCDSR. This information reveals the individual contri- 
butions of the u, d and s quarks to the spin content of A. 
We have found that in the SU(3)i? limit our results for g\ 
are in agreement with expectations based on experiment 
assuming SU(3)i? symmetry while the value we obtain for 
g^ slightly deviates from the empirical one. We have also 
analyzed the isoscalar coupling g'^ with SU(3) f breaking 
effects and have found that the light-quark contributions 
remain mainly unaffected and negatively polarized as in 
the SU(3)i. limit. 



APPENDIX A 

A simple model for the A baryon is that it is con- 
taminated by E baryon together with a tt having orbital 
angular momentum L = 1 in order to conserve parity. In 
the naive quark model the S is composed of a u-d quark 
pair coupled to 5* = 1 and a s-quark with S = 1/2. In 
this picture, the quark-spin configurations of S are given 



as 



|I](+l/2)) = x/273|M](+l),5(-l/2)> 
-yT73|M](0),s(+l/2)), 

|S(-l/2)) = -v/273|M(-l),s(+l/2)) 
+ V^|M](0),s(-l/2)), 



(Al) 
(A2) 



using the appropriate Clebsch-Gordan coefficients. Sim- 
ilarly, the S-TT mixed state is written as 



|E, tt; J = 1/2) = - V273|S(-l/2), L,(+l)) 
+ V^|S(+l/2), L^O)). 



(A3) 



Inserting the spin configurations in Eqs. (|Aip and (|A2[) 
into Eq. (|A3p . we obtain 



(A4) 



IS, tt; J=l/2) = 

2/3|M(-l), s(+l/2), L,(+l)) 
-\/2/3|M(0), s(-l/2), L,(+l)) 

+ V2/3|M(+1), s(-l/2), LAO)) 
-1/3|M](0), s(+l/2),L^(0)). 

It is then straightforward to calculate the spin proba- 
bilities of the u-d quark pair in the S-tt mixed state as 
2/9, 3/9 and 4/9 corresponding to [ud\{+l), [ud\{0) and 
[ud]{—l) configurations, respectively, which results in an 
expectation value of —2/9 for the spin of the u-d pair. 
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